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A low luminosity state in the massive X-ray binary SAX 

J0635+0533 
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Abstract. The X-ray pulsar SAX J0635+0533 was repeatedly observed with the XMM-Newton satellite in 2003-2004. The 
precise localization provided by these observations confirms the association of SAX J0635+0533 with a Be star. The source was 
found, for the first time, in a low intensity state, a factor ~ 30 lower than that seen in all previous observations. The spectrum, 
well fitted by an absorbed power law with photon index ~ 1.7 and Njj = 1.2x 10 22 cm -2 , was compatible with that of the 
high state. The low flux did not allow the detection of the pulsations at 33.8 ms seen BeppoSAX and RXTE data. In view of 
the small luminosity observed in 2003-2004, we reconsider the peculiarities of this source in both the accretion and rotation 
powered scenarios. 
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1. Introduction 

The X-ray source SAX J0635 +0533 was discov ered with 
BeppoSAX in October 1997 dKaaret et al.1 1 1999b during a 
search for counterpart s of unidentified gamma-ray sources 
dThompson et al.lll995l) . Its 2-10 keV flux of 1.2 x 10~ n 
erg cm~ 2 s _1 , hard power law spectrum (photon index ~1.5) 
extending to 40 keV, and positional coincidence with a V = 
12.8 star of Be spectral type, immediately suggested to clas- 
sify SAX J0635+0533 as an accreting high mass X-ray bi- 
nary. The subsequent d iscovery of X-ray pulsations at 33.8 ms 
dCusumano et alj |2000) in the BeppoSAX data, makes this ob- 
ject quite peculiar and raises some problems for the accretion 
scenario. In fact, if the X-ray emission is powered by accretion 

on the neutron star surface, in order to avoid the propeller effect 

I 1 \ — — L 

(see, e.g., [Cjrnrjana_et al. (1995)), the magnetic field in SAX 

J0635+0533 must be about three orders of magnitude smaller 
than expected in a typical high mass X-ray binary. 

The pulse frequencies measured with RXTE in 1999, about 
two years after the source discovery, indicated an orbital mod- 
ulation with a period of about 1 1 days and set a l ower limit on 



the lo ng term spin-down P > 3.8x10 u ss (Ka aret et al 



2000). Such a high P in a rapidly spinning neutron star implies 
a large rotational energy loss, E rot = 10 45 47r 2 P/P 3 = 4x 10 38 
erg s _1 , capable to power the observed X-ray luminosity with- 
out the need of invoking mass accretion. In this interpretation, 
SAX J0635+0533 would resemble other binary systems com- 
posed of a fast pulsar orbiting a massive star, such as PSR 
B 1259-63 dJohnston et all 19921) . in which the X-ray emission 
is thought to originate in the shock between the pulsar's rel- 



ativistic wind and that of the companion star. The failure to 
detec t radio pulsations from SAX J0635+0533 dNicastro et al 



2000), does not rule out this scenario, since beaming and/or ab- 
sorption effects might render the radio emission unobservable. 

Here we report the results of XMM-Newton observations, 
carried out in 2003-2004, during which SAX J0635+0533 was 
detected at a very low flux level, the smallest ever seen from 
this source. 



2. Observations and data analysis 

SAX J0635+0533 was observed by XMM-Newton with ten dif- 
ferent pointings between 2003 September 1 I th and 2004 April 



14 . The three E PIC focal plane cameras dTurner et alj|2001 
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Stru der et alj2001 ) were active during these pointings. The two 
MOS cameras were operated in the standard Full Frame mode 
(time resolution 2.6 s), in order to cover the whole 30' field-of- 
view. The pn camera was operated in Timing mode, providing 
a time resolution of 30 /is, but without imaging information. 
For each observation and focal plane camera the thin filter was 
used. 

For each pointing we retrieved from the XMM-Newton 
archive the pps files produced by pipeline processing system, 
which is operated by the Survey Science Center. Since only 
Timing mode data were obtained with the pn camera, we con- 
centrated on the MOS data to study the source properties. Due 
to the faintness of SAX J0635+0533 during these observations, 
it was not possible to see the 33.8 ms pulsations in the pn cam- 
era data. For each observation, we looked for possible periods 
of high instrumental background, caused by flares of soft pro- 
tons with energies less than a few hundred keV. To this aim, we 
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selected only single and double events (PATTERN < 4) with 
energies greater than 10 keV and recorded in the peripheral 
CCDs. Then, we set a count-rate threshold for good time inter- 
vals (GTI) at 0.5 cts s _1 . By selecting only events within GTIs 
we finally obtained a "clean" event list for each MOS data set. 
The dates and effective exposure times (after soft-proton rejec- 
tion) of the observations are listed in Table Q] 

After merging the event lists of the two MOS cameras, we 
accumulated an image of the field-of-view comprising all the 
observations. This clearly showed a source, with a count rate 
of (4.3 ± 0.3) x 10~ 3 counts s _1 in each camera, at the coordi- 
nates R.A. = 06 h 35 m 18.3 s , Dec. = +05° 33' 06.3" (J2000). 
This position was obtained after the astrometric correction of 
the X-ray coordinates of the detected sources, based on the 
positions o f five optical count erparts found in the Guide Star 
Catalogue dLasker et al.1 120081) . The final uncertainty on the 
source position is 1" at 90% c.l. (including statistical and sys- 
tematic errors). The XMM-Newton localization confir ms the 
association with the star proposed by lKaaret et al. (1999) on the 
basis of an estimated 4% probability of finding by chance a Be 
star in the 30" radius error circle determined with BeppoSAX 
(see Fig. |T}. 

We extracted the source spectra by selecting events in a cir- 
cular region with a small radius (10") in order to minimize the 
background contribution. The background spectra for the two 
MOS cameras were accumulated from large circular areas with 
no sources and radii of 100". We generated ad hoc redistribu- 
tion matrices and ancillary files using the SAS tasks rmf gen 
and arf gen, respectively. In order to ensure the applicabil- 
ity of the x 2 statistics, all spectra were rebinned with a mini- 
mum of 30 counts per bin, and fitted in the energy range 0.3- 
10 keV using XSPEC 11.3.2. We checked that separate fits of 
the two spectra gave consistent results, therefore we analysed 
them simultaneously, in order to increase the count statistics, 
imposing common spectral parameters for the two spectra and 
introducing a free relative normalization factor between them, 
to account for possible differences in the instrument cross- 
calibration. We found a 1.02 cross-normalization factor be- 
tween the MOS2 and the MOS1 camera. 

A good fit was obtained with an absorbed power- 
law (Fig|2|i, yielding a hydrogen column density A^h = 



+0.43 
-0.31- 



(1.17+0 ^) x 10 22 cm" 2 and a photon index Y = 1.74 
The absorbed flux in the energy range 0.2-12 keV is fx ~ 
1.45 x 10~ 13 erg cm~ 2 s -1 , while the corresponding unab- 
solved flux is ~ 2.55 x 10~ 13 erg cm~ 2 s _1 . Acceptable fits 
were also obtained with a thermal bremsstrahlung (kT = 9.2 
keV) and with a blackbody (kT ~ 1 keV). 

To study the source variability, we analyzed separately 
the ten data-sets. For each of them we performed a detailed 
source detection in five energy bands (0.2-0.5, 0.5-1, 1-2, 2- 
4.5 and 4.5-12 keV), applying the same procedure and param- 
eters used by the XMM-Newton SSC to produce the XMM - 
Newton serendipitous source catalogue ( Watson et al. 2009). 
The source detection was performed simultaneously on both 
MOS data sets and using the corresponding exposure maps 
which account for spatial quantum efficiency variations, mirror 
vignetting, and effective field of view. The threshold value for 
the detection likelihood was set to 5 for the likemin parameter 



Fig. 1. Optical image of a 2'x2' field around the position of 
SAX J0635+0533 from the Digitized Sky Survey (B filter). 
The larg e circle is the B eppoSAX error region with a radius 
of 30" (ICusumano et al. 2000). The small circle (2" radius) is 
the XMM-Newton position. North is to the top, East to the left. 




channel energy (keV) 



Fig. 2. Top panel: average spectrum of SAX J0635+0533 with 
the best-fit power-law model. The spectra of the MOS1 and 
MOS2 cameras are shown in black and red, respectively. 
Bottom panel: data-model residuals, in units of a. 



of the SAS task eboxdetect and to 6 for the mlmin parameter 
of the SAS tasks emldetect and esensmap. We emphasize 
that these values imply a rather loose constraint on the source 
significance in order to be detected; they allow to detect even 
very weak sources, with a low (i.e. ~ 3) signal to noise ratio. 
In this way we aim to check if SAX J0635+0533 is even only 
marginally detected in any observation. On the other hand, in 
our detection procedure we used 'ad hoc' energy conversion 
factors (ECF) to convert the measured count rates of the de- 
tected sources (both in each of the five energy bands and in 
the total 0.2-12 keV band) into the corresponding energy flux. 
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They were derived based on the best-fit power-law model of 
the source average spectra. 

Based on this analysis, we found that SAX J0635+0533 
was detected only in six observations, as reported in Table Q] 
The upper limits on the count rate (at the confidence level cor- 
responding to a detection likelihood L = 6) are obtained from 
the sensitivity maps at the source position. 

The long term light curve of SAX J0635+0533 is plotted 
in Fig. [3] In September-October 2003 the source flux varied by 
at least a factor 10. Although the observations are not continu- 
ous, they suggest an outburst lasting about three weeks, with a 
rise time of only a few days to a maximum flux of ~ 5 x 10 -13 
erg cm" 2 s _1 , followed by a similarly rapid decay. A simi- 
lar flux level was observed again six months later. Comparing 
the hardness ratios measured in the different observations, we 
found some evidence for a slight spectral hardening correlated 
with the source intensity (Table [T] Fig.@]i. 
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Fig. 4. Hardness ratio versus source flux. The hardness ratio is 
defined as (H-S)/(H+S), where H and S are the source count 
rates in the hard (H = 2-12 keV) and soft (S = 0.2-2 keV) 
energy ranges. 
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Fig. 3. Light curve of SAX J0635+0533 . The count rates refer 
to the 0.2-12 keV and to the sum of 2 MOS. The data of the first 
two observations have been merged. The upper limits (obtained 
with a threshold in detection likelihood L = 6) correspond to a 
~ 3cr confidence level. 



3. Discussion 

The maximum flux we observed for SAX J0635+0533, 4.2 
x 10~ 13 erg cm~ 2 s _1 in the 2-10 keV range, is a factor > 
30 smaller than that measured at the time of the BeppoSAX 



discov ery in 1997 (1.2 x 10 erg cm 1 s . iKaaret et al 



1999)). To our knowledge, this is the lowest flux ever reported 



for SAX J0635+0533. The upper limits of the XMM-Newton 
pointings of September 2003 imply an even smaller flux. 

The distance of SAX J0635+053 3 is not well constr ained. 
The range 2.5-5 kpc was estimated by Kaaret et al. ( 19991) from 



the properties of the proposed optical counterpart. We also note 
that a distance much in excess 5 kpc is unlikely, considering the 
location of SAX J0635+0533 in the Galactic anti-center direc- 
tion. In the following discussion, where we consider the two 
alternative possibilities for the origin of the observed X-ray 



emission, we conservatively normalize the relevant quantities 
to a distance ds of 5 kpc. 

3. 1. Accretion powered X-ray emission 

Neutron stars accreting from Be star companions constitute the 
large majority of the high mass X-ray binary systems present in 
the Galaxy. Also in view of our improved localization of SAX 
J0635+0533, it is natural to first discuss this possibility. 

Already at the time of its discovery, it was noticed that the 
luminosity of SAX J0635+0533 was relatively small compared 
to classical Be/neutron stars systems. The persistent sources of 
this class have X-ray luminosity of 10 36 — 10 37 erg s _1 . This 
is also the luminosity typically reached during the outbursts 
of transient Be/neutron stars systems, which comprise the ma- 
jority of this population. With the advent of more sensitive 
observations a number of persistent Be binaries with lower 
lumi nosity, ~ 10 34 erg s ~ 1 , have also been discovered (see, 



e.g. iReig & Rochd dl999l) : lLa Palombara & Mereghettil feOOfj. 
2007b ). Our observations indicate for SAX J0635+0533 an av- 
erage luminosity, a few 10 33 d| erg s _1 , much smaller than 
these values, and provide an upper limit as low as 3x 10 32 d| 
erg s _1 in mid September 2003. If we further consider that the 
source could well be closer than 5 kpc, we are faced with an 
even smaller luminosity. 

The new data clearly indicate that SAX J0635+0533 is a 
transient source, but it differs from the other Be systems for 
its low luminosity both during the "high state" and in "quies- 
cence". The non detection in September 2003 allows us to set 
an upper limit on the mass accretion rate of 1.5 xlO 12 g s _1 . 
This limit applies assuming that the accretion flow proceeds 
down to the neutron star surface, which is very unlikely if in- 
deed the neutron star is rotating at 33.8 ms. In the presence 
of the neutron star magnetic field, different scenarios prevent- 
ing, or reducing, the accretion ra te onto the neutron sta r sur- 
face can occur (see, for example, ICampana et al.1 (119981) ). For 
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Table 1. Flux and luminosity values of SAX J0635+0533 in the individual observations. 



Start Observation 


Net exposure 


Absorbed flux (0.2-12 keV) 


Luminosity'"' 


Hardness 


date - UT 


(ks) 


(erg cm 2 s 


) 


(erg s — 1 ) 


Ratio (b) 


2003-09-11 13:22:41 


8.1 


< 6.1 x 10" 


14 


< 3.03 x 10 32 




2003-09-15 18:47:20 


5.1 


< 7.8 x 10" 


14 


< 3.88 x 10 32 




2003-09-18 13:17:46 


6.8 


<r 9 9 v m~ 


13 


< 11.1 x 10 32 




2003-09-21 09-53-21 


3.7 


(1.7 ±0.5) x 10" 


13 


(8 5 ± 5) x 10 32 


0.17 ± 0.18 


2003-09-25 12:45:02 


5.0 


(5.5 ±0.7) x 10" 


13 


(27.6 ± 3.4) x 10 32 


0.48 ± 0.08 


2003-10-08 07:29:55 


2.8 


(3.7 ±0.6) x 10" 


13 


(18.6 ±2.4) x 10 32 


0.29 ±0.13 


2003-10-13 15:30:19 


4.6 


(9.0 ±3.3) x 10" 


14 


(4.5 ± 1.6) x 10 32 


-0.28 ± 0.21 


2003-10-15 22:34:19 


5.9 


(3.2 ± 1.3) x 10" 


14 


(1.6 ±0.7) x 10 32 


-0.06 ±0.15 


2004-03-13 11:44:58 


2.7 


(4.0 ± 0.7) x 10" 


13 


(20.7 ± 3.3) x 10 32 


0.30 ±0.11 


2004-04-14 08:04:48 


4.5 


< 1.6 x 10" 


13 


< 8.0 x 10 32 





In the energy range 0.2-12 keV, corrected for the absorption, and assuming a distance of 5 kpc. 

Ratio between the source count rates in the hard (H = 2-12 keV) and soft (S = 0.2-2 keV) energy ranges, defined as (H-S)/(H+S). 



such a short spin period and low luminosity, the direct accretion 
regime, in which the magnetospheric radius is smaller than the 
corotation radius and hence the magnetic centrifugal barrier is 
open, requires a magnetic field smaller than ~ 10 8 G. This field 
is at least three orders of magnitude lower than that expected in 
a young neutron star with a Be companion. In fact all other ac- 
creting pulsars of this class have much longer spin periods. The 
only exception is the recurrent tran sient A 0538-66 which ro- 
tates at 69 ms ( Skinner et al.ll 1 982b . Note that the pulsations in 
this systems were only detected during a brig ht outburst r each- 
ing a luminosity of ~ 10 39 
plying a magnetic field of ~ 



erg s _1 (ISkinner et al.ll 1982b . im- 
10 11 G. Although smaller than the 
average, such a field is not implausible. Furthermore it is con- 
sistent with the interpretation of the unpulsed quiescent lumi- 
nosity of A 0538-66 (sev eral 10 33 erg s -1 ) as accretion halted 



at the centrifugal barrier (ICampana et al 1 ll995h . 

If SAX J0635+0533 has a typical magnetic field, the low 
luminosity observed with XMM-Newton could be due to mass 
accretion stopped at the magnetospheric radius. Assuming for 
simplicity spherically symmetric accretion, and a neutron star 
with mass 1.4 Mq and radius 10 6 cm, gives in this case an 
X-ray luminosity of - 2 x 10 32 B^ /7 M% 7 ergs s" 1 (Bi 2 
is the magnetic field in units of 10 12 G and Mi 5 the accretio n 



rate in units of 10 15 g s 1 ; see, e.g.. ICampana et al.l (119981) 1 



However, the higher luminosity state observed in the past with 
BeppoSAX and RossiXTE cannot be explained in the same 
way due to the presence of pulsations with a relatively high 
pulsed fraction, which are not expected when the magnetic cen- 
trifugal barrier operates. 

In conclusion, the difficulties already pointed out in inter- 
preting SAX J0635+0533 as a typical accretio n powered Be 
binary (IKaaret et al. 2000; Nica stro et al.l l2000) are reinforced 
by the low luminosity reported here. Of course, the properties 
of this system would fit better in this scenario if the fast peri- 
odicity were disproved by further observations. 

3.2. Rotation powered X-ray emission 

The alternative interpretation of SAX J0635+0533 is that of a 
rotation powered neutron star, whose X-ray emission derives 
from the shock between the relativistic pulsar wind and the 



companion's wind. The large luminosity difference between 
our data and the previous observations could be due to vary- 
ing shock conditions in a very eccentric orbit. However, also 
in this scenario this source would present some peculiar prop- 
erties, compared to the (admittedly few) other systems of this 
kind. The large variations seen in September 2003 are difficult 
to explain if the source was far from periastron, where no big 
changes in the shock properties are expected. Rotation powered 
pulsars in intera cting binaries, such as the alrea dy mentioned 
PSR B 1259-63 dChernvakova et alJl2006ll2009h or the "black 
widow" pulsar PSR B 1957+20 dHuang & Becker! 120071) . have 
X-ray efficiencies i n the range 10~ 4 — 10 2 . On the other hand, 
the P reported by (IKaaret et al ]l2000h . corresponding to a rota- 



tional energy loss E rot larger than a few 10 38 erg s 1 , implies 
a much lower efficiency for SAX J0635+0533. 

4. Conclusions 

Despite SAX J0635+0533 was found by XMM-Newton in a 
very low luminosity state, thanks to the good sensitivity of the 
EPIC instrument, we could derive a precise localization that 
confirms the proposed Be optical counterpart of this source. We 
observed large flux variability when the source was detected in 
September/October 2003 and derived a stringent upper limit of 
3 x 10 32 d 2 , erg s _1 for its luminosity when the source was not 
detected. 

The spectral and flux properties of SAX J0635+0533 are 
consistent with a variable neutron-star binary powered by ac- 
cretion from the Be companion or by the loss of rotational en- 
ergy. Both interpretations imply some peculiarities with respect 
to other known sources. These result from the very short spin 
period and possibly high period derivative reported for SAX 
J0635+0533, which unfortunately we could not confirm owing 
to the source faintness during our observations. 
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